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Our  first  attempts  at  a  quantitation  of  the  mechanism  of 
proton  propagation  in  thin  water  layers  between  phospholipid 
membranes  was  based  on  microsecond  kinetic  measurements  supported 
by  the  steady  state  Debye  Smoluchowski  equation  for  diffusion 
controlled  reactions  as  an  analytic  tool.  These  measurements 
revealed  the  effect  of  the  width  of  the  water  layer  on  the  rate 
of  proton  transfer  reactions,  but  could  not  evaluate  specifically 
the  contribution  of  each  of  the  following  paramaters:  diffusion 
coefficient,  intensity  of  electrostatic  attractions,  stenc 
hindrance  and  uptake  of  proton  by  the  phosphomoiety  of  the 
membrane.  Yet  within  these  restrictions  we  concluded  that  the 
diffusion  coefficient  of  proton  in  thin  water  layer  is  similar, 
or  somehwat  smaller  than  that  in  bulk  water.  We  also  quantitated 
the  rate  constant  of  proton  capture  by  the  phospho  moiety  of 
phosphatidyl-choline  (6xl09JM-1  sec-1)  which  is  '3%  of  the  rate 
measured  with  pyranin  (1) . 

The  lack  of  accurate  expression  of  the  contribution  of 
electrostatic  potential  to  the  rate  of  the  reactio..  forced  us  to 
look  for  a  better  resolving  system. 

The  observation  method  we  selected  was  the  fluorimetric 
monitoring  of  the  geminate  recombination  between  H  and  the 
excited  pyranin  anion  (40  -)  as  it  proceeds  in  the  thin  water 
layer.  The  dynamics  were  followed  with  a  20  ps  resolution  over 
very  short  time  frame  of  2  ns.  The  results  were  analyzed  by  a 
semiquantitative  procedure  which  indicated  that  the  average 
number  of  proton-anion  recombinations  is  a  function  of  the 
physical  characteristics  of  the  reaction  space,  mostly  dominated 
by  the  activity  of  water  (3h20^  laYer  (2)  •  As  activity 
of  water  and  the  average  number  of  recombinations  are 
macroscopic  parameters,  we  looked  for  a  better  system  for  the 
analysis  of  the  results. 

The  present  procedure  satisfies  these  demand^ 
Experimentally  we  follow  the  dynamics  of  the  proton  40*- 
recombination  over^a  longer  time  range  (20  ns)  where  the 
intensity  of  the  40H  fluorescence  decreases  by  almost  4  orders 
of  magnitude,  and  our  analysis  covers  the  whole  dynamic  range. 

The  analysis  of  the  signal  is  achieved  by  N.  Agmon's 
mathematical  procedure  (Department  of  Chemistry,  the  Hebrew 
University  of  Jerusalem;  3,4,5).  The  beauty  and  advantage  of  his 
analytic  formalism  is  the  usage  of  stochastic  arguments  for  the 
description  of  the  events.  These  stochastic  values,  given  by  the 
probabilities  of  events,  are  blended  by  the  numeric  treatment  to 
reproduce  the  time  resolved  dynamics  of  the  observed  population. 
The  stochastic  parameters  are  an  explicit  function  of  physical 
parameters  like  dielectric  constants,  diffusion  coefficients, 
size  of  the  reaction  sphere,  the  width  of  the  aqueous  layer  and 
the  release-capture  of  a  proton  from  all  reactants  present  in  the 
system. 


To  evaluate  the  accuracy  of  the  analysis,  we  tested  the 
whole  procedure,  kinetic  measurements  and  mathematical  analysis, 
on  the  well  defined  system  of  pyranin  dissolved  in  concentrated 
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solutions  of  sucrose.  The  analysis  was  carried  out  using  the 
dielectric  constant  of  the  solution  as  an  adjustable  parameter 
and  was  compared  with  the  values  determined  by  the  classical 
techniques  (6) .  The  results  of  our  calculations  are  depicted  in 
Figure  1.  Thus  the  fast  kinetic  measurements  are  capable  of 
monitoring  the  dielectric  constant  of  the  immediate  vicinity  of 
the  dye. 


FIGURE  1 

The  dependence  of  the  dielectric  constant  of  sucrose 
solution  on  its  molar  concentration.  The  points  shown  in  the 
graph  were  calculated  from  time  resolved  fluorescence  dynamics. 
The  continuous  line  is  drawn  through  the  values  listed  in  the 
International  Critical  Tables  (61. 


The  verification  of  the  validity  of  the  methodology,  as 
tested  in  a  spheric  symmetric  space,  encouraged  us  to  employ  it 
for  the  study  of  the  dynamics  in  the  thin  water  layer  between 
npia  membranes. 


For  the  study  of  proton  transfer 
Agmon's  program  had  to  be  adjusted  to 
factors: 


in  thin  water  layer, 
consider  the  following 


.  ,1,  .Distortion  of  the  spheric  symmetry  of  the  electric  field 

CAarq?f  f£rmed  °n  the  dielectric  boundary  between 
Pld.  (.7J  •  To  account  for  that  unknown  value,  the 
pfJameter.  constant  (eeff)  was  treated  as  an  adjustable 

delavt'  rhe“P°f?,TY  caPture  by  the  phosphomoeties  of  the  membrane 
diJfu^nn1 nE£?f9iatlfn  ??  ^  proton,  leading  to  an  apparent 
coefficient  which  is  smaller  than  that  of  proton  in 

tho  ‘  T°  account  for  that  reaction,  the  binding  of  proton  to 

That  expressed  as  a  stochastic  event  (see  Scheme  I). 

Blall.stlc  treatment  segregated  binding  from  diffusion, 
vilue^  t0  quantitate  the  diffusion  coefficient  by  its  real 


SCHEME  I 


Schematic  presentation  of  the  reaction  space  and  definition 
of  parameters. 

The  pyranin  is  adsorbed  to  the  membrane  and  enclosed  within 
the  reaction  sphere  with  radius  Ro. 

The  proton  ejected  from  the  pyranin  appears  on  the  surface 
of  the  reaction  sphere  at  the  rate  k^,  and  is  readsorbed  with  kr» 
The  propagation  of  proton  between  the  concentric  shells  is  given 
by  the  rate  constants  k.*  and  k*.  .  Tb£'  frequency  element  of  the 

rates  is  controlled  by  the  diffusion  coefficient  and  their  value 
is  function  of  the  gradient  of  the  electrochemical  potential  of 
the  proton  (3) .  The  binding  of  proton  to  the  phosphomoiety  is 
given  by  kas3  and  redissociation  by  k^sg.  Within  the  dotted 
space,  th«  proton  propagates  in  a  3 -aimens ional  space  where 
dv/dr«*f(r  /  beyond  this  range  dv/dr=f{r),  as  in  a  2-dimensional 
space. 


3.  Geometric  constraints  of  the  reaction  space  limit  the 
random  walk  of  proton  to  a  two  dimensional  space  at  the  point 
volnm/nf  h  ^see.  Scheme)  •  In  that  region  the  incremental 
,  ,/U_  v®,rij\s  aS  dv/dr  ~  f  <r>  while  ^  the  range 

incorporated/in  ter^  Wlto 

i0n  °f  “*  probabilities 

n,rf  steric  restriction  of  approach  trajectories  due  to 
iai  masking  of  the  reaction  sphere's  surface  by  contact  with 
and3!3?6;  °n-  the  surface  the  reaction  space 

of f  the  surf  ace  ht  n  actlvity  of  water  (9)  and  the  masking 
ind.nZnJ-n/' t  by  n°n  aqueous  solute.  As  aH?Q  is  almost 
independent  of  sucrose  concentration  we  hacP?o  evaluate 

studies  r^S  °fAh.e  I#asking  of  reaction  sphere.  Our 

almost  HnearP^0t0n  dlssociation  in  sucrose  solutions  yielded  an 

sucrose  in  thIT?  ^  b°tb  *r  a*d  with  the  partial  volume 
analysis  solution  (see  figure  2).  In  accordance,  the 

allowed  dynamics4we  measured  with  the  membranal  system 

a“°led  b°^  *f  and  kr  to  vary  and  the  ratio  of  each,  with 

sSfic ‘hSaSs.r*' w*1"'  u*s  tak,n  »* the  vai“  °£  «» 


FIGURE  2 

The  dependence  of  steric  hindrance  fot  proton  approach  to 
pyranin  anion  on  t]je  partial  volume  of  sucrose  in  the  solution. 
Open  symbols  are  G  values  calculated  from  kf;  solid  symbols  are 
for  values  calculated  from  k,..  1 


Experimental 

Time  resolved  fluorescence  of  pyranin  was  recorded  by  a 
single  photon  counting  at  wavelength  of  iOH  emission  over  a  20ns 
period.  The  number  of  photon  counted  at  the  maximum  is  5000  and 
more  than  1 06  were  counted  in  each  experiment. 

Analysis 

The  signals  were  reconstructed  by  the  program  written  by  N. 
Agmon  and  adapted  to  our  special  conditions  by  Mrs.  E.  Nachliel 
of  our  group.  These  adjustments  in  the  program  were  then  chekced 
and  approved  by  Dr.  Agmon. 

Results 

Figure  3  depicts  the  signal  of  40H*  emission  over  a  20  ns 
period  over  dynamic  range  of  more  than  3  orders  of  magnitude. 

The  reconstructed  dynamics  is  shown  on  both  linear  and 
logarithmic  presentation. 

The  experiments  depicted  in  Figure  3  were  repeated  with  DPPC 
and  DPPC+cholesterol  vesicles  over  a  wide  range  of  osmotic 
pressures.  The  results  were  analyzed  and  the  dependence  of  the 
microscopic  constants  on  the  osmotic  pressure  and  resulting  value 
of  dw  (10)  ate  given  in  Table  I. 
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FIGURE  3 

Fluorescence  decay  dynamics  of  excised  pyranin  (*0H  )  as 
measured  for  the  dye  trapped  between  DPPC  membranes  at  a  distance 
of  9 . 4 A  apart.  The  main  frame  depicts  the  signal  on  a  linear 
scale.  The  insert  presents  the  signal  on  a  logarithmic  Y  scale. 
The  continuous  line  is  the  reconstructed  dynamics,  calculated 
according  to  Agmon,  as  modified  for  a  thin  layer,  using  the 
values  given  in  Table  I. 

Discussion 


Ihe  diffusion  coefficient  of  H+ 


In  our  analysis. 
Vet  it  converged  in  z 
water  (Dh+  =  9.3  x  10- 


[+  was  treated  as  an  adjustable  parameter. 

.  samples  to  the  value  character izina  bulk 
cirrysec) .  ' 


bfs's  of  thG  many  experiments  we  have  carried  out,  we 
reject  the  notion  of  Prats  et  al.  (11)  that  the  ordering  of  water 

the  d%fUr-  ^°f  nielnbran®  U2)  has  any  accelerating  effect  on 
the  diffusion  of  proton  (ll). 

Dielectric  constant 


Hiff„i'e»ef£eCt^e  dielectric  constant  calculated  for  the 
di  fusion  space  is  significantly  lower  than  that  of  bulk  water 

dielectr?  ^  <Tontribution  multiple  image  charges  on  the 

the1^^  imniI?darttS.-0f  thS  tW°  menbranes  (?)•  Examination  of 
:?*  ™luf  implies  that  upon  compression  the  intensivity  of  the 

d??f09e'1"'ari.Charges  is  increased .  We  also  note  a  marked 
latter e/*Ce  betwee"  DPPC+cholesterol  surfaces.  In  the 

ZZZ  €eff  ls  smaller,  indicating  a  better  adsorption  of  the  dve 
to  the  surface  (see  also  Gould,  13) .  y 

M.  ^  ®gmbrane:  The  probability  of  proton  capture 
f  phosphomoiety  on  DPPC  membrane  is  '2%  with  respect  ,o  the 
rate  of  recapture  by  pyranin  anion.  This  ratio  is  in  agreement 
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FIGURE  4 

The  effect  of  the  geometry  oJ_reactlon  space  on  the  dynamics 
of  the  reaction  between  H+  and  to 


The  curves  depict  the  expected  fluorescence  decay  of  »01i 
proceeding  in  a  water  layer  between  DPPC-cholesterol  membranes 
(under  no  external  pressure)  as  the  width  of  the ,®5fu®ous eimnet 
(d„)  is  decreased  from  30A  (lowermost  curve)  to  2 OA  (uppermost 
curve) .  As  seen  in  the  figure,  even  the  first  decrament  in  d  (by 
3£A)  already  shifts  the  curve  out  of  the  range  set  by  the 
experimental  noise. 


with  our  previous  analysis,  based  on  Debye  Smoluchowski  equation 
(1)  which  set  the  ratio  as  "3%. 

The  reactivity  of  surface  groups  increases,  as  expected, 
with  the  constriction  of  the  space. 

The  reactivity  of  DPPC  moieties  increases  when  separated  by 
cholesterol  molecules.  This  may  be  interpreted  as  a  reflection  of 
the  larger  separation  between  phospho  moieties  (14)  which 
increase  their  accessibility  to  H+. 

Geometry  of  the  reaction  space ■  The  geometry  of  the  space  is  a 
built-in  feature  in  our  computations.  A  demonstration  of  its 
effect  on  the  dynamics,  insulated  from  the  other  variables,  is 
attainable  by  simulated  dynamics  (see  Figure  4).  This  figure 
demonstrates  that  even  a  2 . 5A  change  in  dw  already  distorts  the 
dynamics  beyond  the  margins  set  by  the  experimental  noise. 

It  should  be  stressed  that  the  width  of  the  water  layer, 
derived  from  X-ray  studies  of  multllamellar  vesicles,  is  still 
under  debate  (10,14).  The  estimation  derived  from  the  kinetics  of 
reaction  within  the  aqueous  layer  may  provide  an  independent 
procedure  to  evaluate  the  physical  characteristics  of  this 
ultrathin  layer. 


Steric  restriction;  The  rate  of  proton  production  on  the  surface 
of  the  reaction  sphere,  and  the  rate  at  which  it  is  adsorbed 
there,  were  introduced  as  independent  variables  in  the  program. 
As  seen  in  Table  I,  both  values  are  similarly  affected  by  the 
osmotic  constriction  of  the  aqueous  layer.  We  interpret  this  as 
enhanced  "adhesion"  of  the  pyranin  molecule  to  the  lipid  membrane 
and  masking  the  surface  feo/mthe  aqueous  matrix. 

The  lower  value  of  <T  recorded  for  pyranin  between 
DPPC+cholesterol  membrane  is  thus  in  accord  with  the  lower  value 
of  eeff  measured  with  that  system. 


TABLE  I 

The  effect  of  osmotic  pressure  on  the  dynamics  of  geminate 
recombination  in  water  layer  between  phospholipid  membranes 

Sample  osmotic  dw^1^  kf  kr  kas  €eff 

^dyn/cm^)  (A)  (given  in  109  x  S-1) 


Bulk  water  NA 

NA 

7 

5.5 

NA 

1 

80 

DPPC  vesicles 

in  water  O 

19.6 

3.6 

2.8 

0.15 

0.48 

50 

0.4M  sucrose  l.lxio7 

12.8 

2.7 

2.2 

0.2 

0.4 

50 

0.6M  sucrose  1.8X107 

11.8 

2.7 

2.1 

0.2 

0.4 

50 

1.0M  sucrose  3.57xl07 

10.5 

2.7 

2.1 

0.22 

0.4 

50 

1.5M  sucrose  6. 3 5x10 7 

9.4 

1.2 

1.0 

0.22 

0.18 

50 

1.75M  sucrose8.8  xlO7 

8.8 

1.15 

0.9 

0.17 

0.16 

52 

DPPC+cholesterol 
in  water  0 

28 

5 

4 

0.15 

0.71 

54 

0.5M  sucrose  1.47xl07 

15.8 

1.8 

1.2 

0.35 

0.23 

40 

1.0M  sucrose  3.57xl07 
1.5M  sucrose  8.8xl07 

12 

1.4 

1.1 

0.38 

0.20 

42 

10.2 

1.0 

0.8 

0.32 

0.14 

42 

Data  of  Parsegian 

Taken  as  ratio  of 

(10)  and 

Ms,s 

refs. 

therein 

or 

kf/V 

*4>; 
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